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ABSTRACT 

Emission of two short hard X-ray bursts on 2009 June 5 disclosed the existence of a new soft 
gamma-ray repeater, now catalogued as SGR0418H-5729. After a few days. X-ray pulsations 
at a period of 9.1 s were discovered in its persistent emission. SGR0418H-5729 was monitored 
almost since its discovery with the Rossi X-ray Timing Explorer (2-10 keV energy range) and 
observed many times with Swift (0.2-10 keV). The source persistent X-ray emission faded 
by a factor ^10 in about 160 days, with a steepening in the decay about 19 days after the 
activation. The X-ray spectrum is well described by a simple absorbed blackbody, with a 
temperature decreasing in time. A phase-coherent timing solution over the ^160 day time 
span yielded no evidence for any significant evolution of the spin period, implying a 3a upper 
limit of 1.1 X 10^^"^ s s^^ on the period derivative and of ^3 x 10^'' G on the surface dipole 
magnetic field. Phase-resolved spectroscopy provided evidence for a significant variation of 
the spectrum as a function of the stellar rotation, pointing to the presence of two emitting 
caps, one of which became hotter during the outburst. Finally, a deep observation of the field 
of SGR0418H-5729 with the new Gran Telescopio Canarias 10.4-m telescope allowed us to 
set an upper limit on the source optical flux of i' > 25.1 mag, corresponding to an X-ray-to- 
optical flux ratio exceeding 10'', consistent with the characteristics of other magnetars. 

Key words: pulsars: general - stars: neutron - X-rays: individual: SGR 0418H-5729. 



1 INTRODUCTION 

Soft gamma-ray repeaters (SGRs) form a small class of sources 
which emit short (~0. 1 s) bursts of soft gamma-rays with peak lu- 
minosities of ^10*^-10*^ erg s~^. Originally believed to be a sub- 
class of the gamma-ray bursts (GRBs), they were soon recognised 
as a different population because their bursts repeat and are quite 
soft, at variance with canonical GRBs. Apart from the flaring activ- 
ity, sometimes culminating in dramatic, very energetic giant flares 
with peak luminosity exceeding 10*'' erg s^^ (as it was the case of 



the 2004 event from SGR 1 806-20: iHiTrlev et alj|2005h . SGRs dis- 
play persistent emission from 0.1 to hundreds of keV, with pulsa- 
tions of several seconds (P ~ 2-9 s), rapid spin down (P ~ 10~^^ 
s s^^), and with significant variability in period derivative, pulse 
shape, flux, and spectral shape. 

This phenomenology has been interprete d within the con- 
text of the magnetar m odel (Paczvnski 1992; Du ncan & ThompsonI 
ll992l ; lThomDson & Du ncan 1995, 1996). SGRs are in fact believed 
to be young, isolated neutron stars endowed with an ultrahigh mag- 
netic field {B ~ lO^^-lO^^ G), which is thought to be the ultimate 
energy reservoir for their activity. It is now commonly accepted that 
anomalous X-ray pulsars (AXPs), a peculiar class of pulsating X- 
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ray sources, are also magnetars, in view of the many similarities 
in the ir p ersistent and burs ting emission (see I Woods & ThompsonI 
I2OO6I andlMereehett ill2008l for recent reviews)0 

On 2009 June 5, a new SGR, namely SGR 04 18+5729, was 
discovered thanks to the detection of two short X-ray bursts with 
the Fermi G amma-ray Burst Monitor (GBM) and other hard X-ra y 
instruments jvan der Horst et ai] |2009l ; Ivan der Horst et al.ll2010h . 
After a few days, on 2009 June 10, X-ray pulsations at a period 
of 9.1 s were discovered in the persistent emission with Rossi X- 
ray Timing Explor er (RXTE), further supporting the ma gnetar na- 
ture of the source iGQgus. Woods & Kouveliotoull2009h . Ground- 
based follow-up observations in optical and infrared wavebands 
(discussed in the following), as wel l as at radio wavelen gths (with 
the 100-m Green Bank Telescope; iLorimer et al.ll20()9h . failed to 
detect the source. 

Here we report on the analysis of a series of observations car- 
ried out in the soft X-ray range (1-10 keV) by the X-Ray Telescope 
aboard Swift and by the Proportional Counter Array instrument of 
the RXTE satellite, collected in the period from 2009 June 1 1 to 
2009 November 24. We also present the results of a very deep op- 
tical observation performed on 2009 September 15 with the new 
Gran Telescopio Canarias 10.4-m telescope on La Palma (Canary 
Islands, Spain), which is presently the largest single-aperture opti- 
cal telescope in the world. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 Swift 

The X-Ray Telescope (XRT: [Burrows et^l2005h on-board Swift 
uses a CCD detector sensitive to photons between 0.2 and 10 keV 
with an effective area of about 110 cm^ and a field of view (fov) 
of 23' diameter. Several observations of SGR 0418+5729 were per- 
formed following its discovery, in both photon counting (PC) and 
windowed timing (WT) modes. In PC mode the entire CCD is read 
every 2.507 s, while in WT mode only the central 200 columns are 
read and only one-dimensional imaging is preserved, achieving a 
time resolution of 1.766 ms (see . Hill et al. 2004 for more details). 
Table [T] (see also Fig. [TJ reports the log of the observations that 
were used for this work. 

The data were processed with XRTPIPELINE (version 0.12.3, 
in the HEASOFT software packag^l 6.6), and standard filtering 
and screening criteria were applied by using FTOOLS tasks.We ex- 
tracted the PC source events from a circle with a radius of 20 pix- 
els (one pixel corresponds to about 2'.' 36) and the WT data from 
a 40 X 40 pixels box along the image strip. To estimate the back- 
ground, we extracted PC and WT events from source-free regions 
far from the position of SGR 0418+5729. 

For the timing analysis, the photon arrival times were con- 
verted to the Solar system barycentre with the BARYCORR task 
and using the Chandra source position (RA = 04'' 18"" 33= 867, 
Dec. = +57°32'2 2'.'91, epoch J2000; 0'.'35 accurate) provided by 
I Woods et alj ( l2009h . For the spectral fitting, the ancillary response 
files (arf) were generated with XRTMKARF, and they account for 
different extraction regions, vignetting and point-spread function 
corrections. We used the latest available spectral redistribution ma- 
trix (rmf) in CALDB (vOll). 

^ Ten AXPs and six SGRs are confirmed, and there are a few candidates; 
see catalog at http://www.physics.mcgill.ca/~pulsar/magnetar/main.html. 
^ See |http://heasarc. gsfc.nasa.gov/doc s/software/ftools/ftools-menu.html| 



Table 1. Journal of the 5n7/f/XRT observations of SGR 0418+5729. The 
complete observation sequence number is prefixed by 00031422 followed 
by the three digit segment number given here (e.g. 00031422001). 



Sequence/mode 


Start/End time (UT) 


Exposure" 




yyyy-mm- 


-dd hh:mm 


(ks) 


001 / PC 


2009-07-08 20:52 


2009-07-08 22:54 


2.9 


002 / PC 


2009-07-09 00:07 


2009-07-09 19:29 


10.6 


003 / PC 


2009-07-10 00:17 


2009-07-10 07:01 


5.6 


004 / WT 


2009-07-12 00:29 


2009-07-12 15:04 


7.1 


006 / WT 


2009-07-15 15:16 


2009-07-15 23:40 


7.7 


007 / WT 


2009-07-16 00:55 


2009-07-16 23:36 


16.4 


008 / PC 


2009-09-20 21:11 


2009-09-21 23:04 


9.4 


009 / PC 


2009-09-22 00:46 


2009-09-22 23:28 


7.6 


010 /PC 


2009-11-08 00:38 


2009-11-08 23:03 


15.1 



" The exposure time is usually spread over several snapshots (single con- 
tinuous pointings of the target) during each observation. 

2.2 RXTE 

The Proportional Counter Array (PCA; llahoda et all 11996 *) on- 
board RXTE consists of an array of five coUimated xenon/methane 
multianode Proportional Counter Units (PCUs) operating in the 2- 
60 keV energy range, with a total effective area of approximately 
6500 cm^ and a full width at half maximum (fwhm) fov of about 
1°. The 46 RXTE/PCA observations of SGR 0418+5729 (obs. ID: 
94048) reported here span from 2009 June 1 1 to 2009 November 
24 (see Fig.[Tll. The exposures range from 0.6 ks to 13.6 ks, for a 
total of 194.2 ks. 

Raw data were reduced using the FTOOLS package. To study 
the timing properties of SGR 04 18+5729, we restricted our analy- 
sis to the data in Good Xenon mode, with a time resolution of 1 /is 
and 256 energy bins. After correction to the Solar system barycen- 
tre with the FXBARY task, the event-mode data were extracted in 
the 2-10 keV energy range from all active detectors (in a given 
observation) and all layers, and binned into light curves of 0.1 s 
resolution. 

For the spectral analysis, we have only made use of the Stan- 
dard 2 data (time resolution of 16 s and 128 energy bins) from 
PCU-2, which is the best-calibrated proportional counter unit in the 
PCA, and spectra were extracted including only top-layer events. 
Good time intervals were selected to ensure stable pointing and fil- 
ter out electron contamination events, following the standard pro- 
cedure for faint sources. Response matrices (combining rmf and 
arf files) were generated using PCARSP. Background subtraction 
for PCA data is performed using a synthetic model that accounts 
for the time-variable particle background and for the average cos- 
mic X-ray background (see Jahoda et al. 2006 for details). Possi- 
ble Galactic diffuse emission remains in the background- subtracted 
data, as well as emission from any X-ray source that occurs in the 
instrument fov. 



3 SPECTRAL ANALYSIS 

Because of its nonimaging nature and wide field of view, the 
RXTEfPCA instrument is not well suited for the spectral study 
of relatively faint and soft sources like SGR 0418+5729. The 
5vvi/f/XRT instrument, instead, offers a good sensitivity and ensures 
proper background subtraction (also in the case of WT observa- 
tions). For this reason we started the spectral analysis by working 
only with the Swift/XKF data. 
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Table 2. Swift spectral results for SGR 04 18+5729 for the single blackbody model. Errors are 
at Icr confidence level for a single parameter of interest. The A^h value is l.liE5'2 x 10^^ 



-2 



Sequence 


Time since to 


kT 


Radius'" 


Absorbed flux'^ 




(days) 


(keV) 


(km) 


(10~" erg cm-2 s'^) 


001-003 


33.004-34.430 


0.96 ±0.01 


0.62 ±0.02 


1.31 ± 0.03 


004-007 


36.156^1.121 


0.94 ±0.01 


0.65 ±0.02 


1.35 ±0.02 


008-009 


107.020-109.113 


0.84 ±0.02 


0.43 ±0.02 


0.37 ± 0.01 


010 


155.164-156.098 


0.83 ±0.03 


0.33 ±0.02 


0.21 ± 0.01 



" We assumed as to MJD 54987.862, which is the time of the detection of the first burst. 
^ The blackbody radius is calculated at infinity and for a distance of 5 kpc. 
In the 1-10 keV energy range. 



3.1 Swift/XRT 

The observations carried out within a few days and in which the 
count rates were similar were added in order to accumulate suffi- 
cient statistics for meaningful spectral fits; as a result we obtained 
four spectra: 001-003, 004-007 (WT data), 008-009, and 010 (see 
Table [TJ. The spectral channels were then grouped to have bins 
with a minimum number of 30 photons. The spectral fi tting, with 
the XSPEC analysis package (version 12.4: 1 Arnaudll99d) . was per- 
formed in the 0.7-10 keV energy range, since the very few counts 
from SGR 0418+5729 made the data below 0.7 keV useless. 

We first tried fitting simple models to the spectra: a power law, 
a blackbody, a double blackbody, and a blackbody plus power law, 
all modi fied for inte rstellar absorption. The abundances used are 
those of lAnders & G revesse (1989) and photoelectric abs orptio n 
cross-sections are from Balucinska-Church & McCammonI jl992h . 
We fit all spectra simultaneously with the hydrogen column den- 
sity tied up between the four datasets. While the power-law model 
is disfavoured (xt = 2.10 for 384 degrees of freedom [dof]), 
the single blackbody provides an acceptable fit (xt = 1.16 for 
384 dof); see Table |2] for the best-fitting parameters. Both double 
blackbody and blackbody plus power law yield similarly good fits 
ixl values of 1.10 for 376 dof and 1.12 for 376 dof, respectively), 
but the parameters of the additional spectral components are rather 
unconstrained. Therefore we prefer the simplest model, a single 
blackbody with temperature kT changing from ~1 keV to ~0.8 
keV (the derived absorbing column is A^h — 1-1 x 10^^ cm^^). 
The corresponding observed flux decreases from ~1.3 x 10~^^ 
erg cm"^ s~^ to ~2 x 10"^^ erg cm"^ s ~^ (1-10 keV). These re - 
sult s are consistent with those reported bv lCummings et alj J2009l) 
and lGogiis. Kouveliotou & Woods! ( l2009h , which were obtained us- 
ing a smaller dataset. 

Although the present spectra do not allow us to discriminate 
among more complex (e.g. two component or more physical) spec- 
tral models, magnetar spectra in the ~1-10 keV range are typi- 
cally described by a combination of a soft, thermal component (a 
blackbody) and a harder one. The latter is often represented by a 
non-thermal power-law tail, although in some sources a second, 
hotter bl ackbody provides an equally good (or better) interpreta- 
tion ( e.g. lGotthelf & Halpeml200a ; lTiengo. Esposito & Mereghettil 
I2OO8I) . In the magnetar model, the power-law tail naturally arises 
from resonant cyclotron scattering (RCS) of seed thermal pho- 
tons on magnetospheric currents ( Thompson, Lyutikov & Kulkami 
l2002h . In the light of this, we attempted a simultaneous fit of all 
spectra w ith two different RCS mod els: the o ne-dimensional (ID) 
model bv lLvutikov & Gavriill j2006l) (see also lRea et J]|2008h and 
the NTZ model, based on three-dimensional Monte Carlo simula- 



tions dNobili, TuroUa & Zandl2008l ; IZane et alj2009l) . 

We tried several fits, first leaving all parameters free (except 
for the column density), then assuming that either the surface tem- 
perature or some of the magnetospheric parameters are the same in 
all observations (the optical depth in the ID model and the twist 
angle in the NTZ model). In all cases, although some fluctuations 
are allowed in the magnetospheric properties, we found that the 
best fits (xt = 1-15 for both models) are obtained by invoking 
scenarios in which the major change among the different observa- 
tions is in terms of a decrease in the temperature of the thermal 
seed photons (from ^--^0.9 keV to ~0.75 keV in the ID model and 
from ~1 keV to ~0.9 keV in the NTZ fits) accompanied by a de- 
crease in the emitting surface area (of a factor ~4 and ^--^7 in RCS 
and NTZ, respectively). The reprocessing of the thermal radiation 
by the magnetosphere is always quite modest and these results are 
consistent with the single blackbody fits discussed above. Overall, 
this favours a scenario in which heat is deposited in a spot on the 
neutron-star surface, following the bursting activity, that cools off 
and shrinks as time elapses. 



3.2 RXTE/PCA 

The analysis of the PCA spectra was mainly aimed at having flux 
measurements over the entire data span. We thus adopted the ab- 
sorbed blackbody model and fitted it to all the PCA spectra to- 
gether; the parameters were left free to vary, except for the absorp- 
tion column density that was fixed at the value measured with Swift. 
Spectral channels having ene rgies below 3 keV (standard for RXTE 
spectral analysis - see, e.g., IPark et aill2004h and above 10 keV, 
owing to the very low signal-to-noise ratio, were ignored. This re- 
sulted in an acceptable fit (xi ~ 0-78 for 644 dof) with spectral pa- 
rameters similar to those reported in Table [2] (see the bottom panel 
of Fig.[T] but see also below for indication of a likely background 
contamination in these data). Again, there is a clear indication of 
an overall spectral softening during the decay. 

We plot the resulting long-term light curve in the top panel 
of Fig. [T] (the last three observations provided only upper limits). 
The flux decays as a broken power law. Fixing t = al the time 
of the detection of the first bii rst (2009 June 5 at 20:40:48.883 
UTC; Ivan der Horst et al.ll2010l) , the break occurs at 17.7 ± 0.7 
days, when the index changes from qi = —0.27 ± 0.02 to 
Q2 = —0.93 ± 0.04 (here and in the following all uncertainties are 
at 1(7 confidence level). It is important to note that the fit param- 
eters depend on the assumed time origin, as well as on a possible 
steady component of the SGR flux. The uncertainties introduced by 
making different (reasonable) assumptions are larger than the sta- 
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Days past MJD 54987.862 

Figure 1. Top panel: the outburst decay of the persistent X-ray flux of 
SGR 04 18+5729 as seen by RXTE (black filled triangles) and Swift (blue 
filled circles). We assumed as t = MJD 54987.862, which is the time of 
the detection of the first burst. The fluxes are absorbed and in the 2-10 keV 
energy range; the down-arrows indicate upper Hmits at 3(t confidence level. 
The broken power law describing the RXTE flux decay is superimposed in 
red (solid line); we also plotted (red, dashed line) the simultaneous fit of 
the Swift and RXTE data after subtracting a constant flux from the latter 
(see Section [33t . Bottom panel; evolution of the chai'acteristic temperature 
of the blackbody model inferred from the spectral fitting. The power-law 
model describing the decay of the blackbody temperature (see Section[8j is 
plotted in red. The vertical dashed fines indicate the epoch of the observa- 
tion with the Gran Telescopio Canarias (see Section|6). 



tistical errors reported here by a factor ~5. 

A comparison with the Swift flux measurements indicates a 
substantial discrepancy between the values derived from the XRT 
and the PCA instruments. This is apparent in Fig. [T] where we 
plotted the results from the individual SwiftfXKI observations. 
While the uncertainty in the relative calibration of the two in- 
struments (which is 10-20% in the energy band considered here; 
iKirs ch et al. 2005) cannot account for the mismatch, the fact that 
the difference between the Sw//i'/XRT fluxes and the best-fitting 
curve of the RXTE data is consistent with being constant (Ai^ = 
(1.7±0.1) X 10~^^ erg cm~^ s~^) suggests the presence of a back- 
ground contribution to the RXTE data. After subtracting AF from 
the RXTE fluxes, we fit a broken power law to the combined RXTE 
and Swift data, obtaining slightly different parameters for the de- 
cay: a'l = -0.28±0.02, a'2 = -1.17±0.02, and 19.2±0.4 days 
for the break epoch. 

With respect to the flux discrepancy, we note that it is consis- 
tent with the fluctuations of the cosmic X-ray background, which 
in the PCA produce eff ects of the order of 4 x 10~^^ erg cm~^ s~^ 
in the 2-10 keV band jjahoda et ej]|2006t) . We tried to better as- 
sess the situation by examining the X-ray images of the field. 
Apart from SGR 04 18-1-5729, eight point sources are detected with 
a signal-to-noise above 3 in the 2-10 keV band within the 23'- 
diameter fov of the 5w;yi'/XRT detector. Assuming for the absorp- 
tion the total Galactic value in the direction of SGR 04 18-1-5729 
(~5.6 X 10^^ cm^^, see Section [§}, we estimate that their total 
contribution to the flux observed in RXTE for a blackbody spec- 
trum with kT = 0.2 keV is Ri(1.5-2.5) x 10"" erg cm" ^ s"\ 
~(4-6) X 10^^'^ erg cm"^ s"^ for a power-law spectrum with pho- 
ton index F = 2. Four sources outside the XRT fov, but still in the 
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Figure 2. Top panel: 2-10 keV RXTE/PCA (binned at 10 days) and 
SwiftfXKT pulse phase evolut i on wi th respect to the period measured by 
iGogUs. Woods & Kouveliotoij j2009h as a function of time. The red solid 
line represents the phase-coherent timing solution (discussed in Section|4j. 
Bottom panel: time residuals with respect to the timing solution. 



PCA on e, are present in th e ROSAT All-Sky Survey Faint Source 
Catalog IS ges et al ] |2000l : see also SectionO; using the same as- 
sumptions as above, their 2-10 keV flux is estimated to be of ~(3- 
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erg cm s in the case of the blackbody spectrum. 



!(2-5) X 10 erg cm s for the power-law spectrum. 



4 TIMING ANALYSIS 

Both RXTE/PCA and Swift/XRT data were used to study the tim- 
ing properties of SGR 0418-1-5729. We checked all observations for 
the presence of bursts by a careful inspection of the light curves 
binned with different time resolutions, but none was found. We then 
started the analysis of the coherent pulsations by inferring an accu- 
rate period measurement. In order to do this, we folded in the 2-10 
keV band t he PCA (binned at 10 days) and XRT data at the period 
reported bv lGogiis, Woods & Kouveliotoul ( 12009 ) and studied the 
phase evolution within the observation by means of a pha se-fitting 
techn ique (details on this technique are given in DaU'Osso et al.l 
l2003h . Given the intrinsic variability of the pulse shape as a func- 
tion of time, we did not make use of a pulse template to cross- 
correlate with. We inferred the phase of the modulation by fitting 
the average pulse shape of each observation with a number of si- 
nusoids, the exact number of which is variable and determined by 
requesting that the addition of a further (higher) harmonic is not 
statistically significant (by means of an F-test) with respect to the 
null hypothesis. 

The result of this analysis is reported in Fig. (2] The fit of the 
pulse phases with a linear component results in a reduced xt = 1-2 
for 12 dof. The resulting best-fit period is P = 9.078 388 4(2) s 
(epoch 54993.0 MJD). This timing solution implies a root mean 
square (rms) variability of only 0.03 s. The inclusion of a quadratic 
term is not statistically significant and we derive a 3a upper limit 
on the first period derivative of |P| < 1.1 x 10"^^ s s"^. 

The light curve of SGR 04 18-1-5729 in the 2-10 keV band 
folded at our best period is shown in Fig. [3] The pulse morphol- 
ogy is complex, showing two asymmetric peaks per cycle, the sec- 
ond of which exhibits two subpeaks. Because of the problematic 
background subtraction in the RXTEfPCA event-mode data and the 
possible contamination discussed in Section [T2l here we will avoid 
computing pulsed fractions and hardness ratio, which are meaning- 
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Figure 3. 64-bin RXTE pulse profiles before (left panel) and after (right 
panel) the break (see Fig.[T]and Section IX2t . The colours correspond to dif- 
ferent energy bands: total (2-10 keV, top, in black), soft {2-A keV, middle, 
in red), and hard (4—10 keV, bottom, in blue). 

ful only for background-subtracted data. It is nevertheless evident 
that the pulse shape changes as a function of energy; one of the sub- 
peaks is in fact prominent in the hard (4—10 keV) band and nearly 
ab sent in the soft {2-A keV) b and. Furthermore, as already observed 
bv iKuiper & HermsenI ( l2009h . the pulse profile evolves with time. 
This is apparent in Fig. [3] where we show the profiles over the ini- 
tial phase of the outburst and after the break at about 19 days in the 
decay curve (see Fig.[T]and Section [J!2l ). 

SwiftlXKI observations have lower statistics and shorter time 
coverage, and in PC mode the frame time (2.5 s) does not allow to 
resolve the pulse fine structures. However, the XRT data have the 
advantage of providing safe background subtraction and thus the 
measure of the pulsed fractions. To increase the statistics, we com- 
bined the observations as in Section |3] The rms pulsed fraction^ 
in the 0.7-10 keV band in the four segments in chronological order 
are (23 ± 2)%, (34 ± 1)%, (26 ± 4)%, and (39 ± 6)%. Apart from 
an overall increasing trend of the pulsed fraction, which is only 
marginally significant (i.e. <3(t), the most substantial difference is 
between the values of the first two segments. We caution, however, 
that the poorer detail resolution of the PC data, with respect to the 
WT data of the segment 004—007, can result in a smearing of the 
pulse profile (and hence in a reduction in the measured pulse frac- 
tion) making it difficult to compare the values. 

The Swift folded light curves confirm the variability of the 
pulse profile with energy; in particular, this is clearly visible in the 
high resolution WT data shown in Fig. |4] where the pulsed frac- 
tion is (42 ± 4)% in the soft (1.5-2.5 keV) band and (29 ± 3)% 
in the hard (2.5-10 keV) band. To assess the statistical significance 
of the pulse shape variation, we com pared the data us ing a two- 
sided Kolmogorov-Smimov test (see IPress et al.l 0*9921 and refer- 
ences therein). The result shows that the difference between the 
soft and hard energy ranges is highly significant: the probability 



^ We define the pulsed fraction as 4 ST- 1 ~ V)'^\ • where n is 
the number of phase bins per cycle, yi is the number of counts in the i-th 
phase bin, and y is the mean number of counts in the cycle. 
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Figure 4. 32-bin SwiftlXKV pulse profiles (background-subtracted) of 
SGR 0418+5729 in the 1.5-2.5 keV (red) and 2.5-10 keV (blue) energy 
ranges. The vertical dashed lines indicate the phase intervals used for the 
phase-resolved spectroscopy (see Section|5). 



that the two profiles come from the same underlying distribution is 
in fact ~1. 2 x 10"". 



5 PHASE-RESOLVED SPECTROSCOPY 

To study in detail the spectral variations of SGR 0418+5729 with 
its pulsation phase, we extracted six spectra from the SwiftlXKI 
data in WT mode (obs. 004-007) according to the phase intervals 
indicated in Fig. |4] We could not use the Swiftl'KSH data in PC 
mode because the time resolution of this operating mode is larger 
than the width of the chosen phase intervals. On the other hand, the 
RXTEfPCA data are not suited either for this kind of analysis due to 
the large systematic uncertainties affecting the background subtrac- 
tion and calibration of the PCA data with sufficient time resolution. 

A simultaneous fit of these spectra with an absorbed black- 
body model with only the normalisation allowed to vary gives an 
acceptable fit (x^ ~ 1.15 for 365 dof), with spectral parame- 
ters consistent with the phase-averaged spectrum. Allowing also 
the blackbody temperature to vary, a slightly better fit is obtained 
(xt = 1.09 for 360 dof), with the highest temperatures in phase 
D and F (fcT ^ 1.04 keV) and the lowest one during phase C 
{kT ~ 0.87 keV). 

We then concentrated on the two phase intervals with the 
largest differences, that are also apparent in the RXTE and Swift 
pulse profiles (Figures [3] and |4]l: phases C and D. A simultane- 
ous fit of the two spectra with a single blackbody with the same 
temperature gives xt ~ 1-26 for 96 dof, with clearly structured 
residuals (see Fig. [5)- A slight improvement in the fit is obtained 
allowing for a different blackbody temperature in the two spectra 
ixl = 1.12 for 95 dof), that are kT = 0.90l^;J^3 keV for phase 
C and kT = imto'ol keV for phase D. The addition of a second 
blackbody gives an even better fit (xt = 0.97 for 91 dof), with the 
following best-fit parameters (TVh = (4 ± 1) x 10^^ cm~^ for both 
spectra and assuming a distance of 5 kpc): kTi = 0.31+o og keV, 
Ri = 3.4lJ;^ km, fcTa = 0.90 ± 0.02 keV, R2 = 0.791^;!^^ 
km for phase C and kTi = 0.23 ± 0.02 keV, Ri = 6 ± 2 
km, kT2 = l-14lo:o4 keV, R2 = 0.43 ± 0.04 km for phase 
D. A similarly good fit (xl = 0.99 for 94 dof) but with a re- 
duced number of free parameters is obtained by simultaneously 
fitting the two spectra with the same single absorbed blackbody 
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Figure 5. Phase-resolved spectroscopy of the SwiftDiKY observations in 
WT mode. Top panel: data from the phase intervals C (see Fig.|4), in black, 
and D, in red. We also plotted (solid lines) the best-fitting blackbody models 
(with the temperatures forced to a common value, see Section |5). Bottom 
panel: residuals from the models in units of standard deviations. 

model (iVn = (7 ± 5) x 10^° cm~^ kT = 0.88 ± 0.03 keV 
and R — 0.83 ± 0.05 km) and a broad Gaussian absorption line 
(E = 2.0 ± 0.1 keV, as = 1.0 ± 0.1 keV and equivalent width of 
1.1 keV) only in the spectrum corresponding to phase D. 
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Figure 6. 20 x 15 arcsec^ section of the i'-band image of the field of 
SGR0418-H5729 as observed by the OSIRIS camera of the GTC. North is 
up, east to the left. The radius (O'.'S) of the error circle of SGR0418-H5729 
(in red) accounts for the 0'.'35 uncertainty on the Chandra position and for 
the accuracy of our astrometry calibration (see Section|6). Five objects are 
detected in the field (labeled with capital letters); the observed apparent 
magnitudes of the brightest sources are: i' = 19.536 ± 0.004 (A), i' = 
20.117 ± 0.007 (B), and i' = 20.290 ± 0.007 (C). The two sources (B 
and D) visible northeast of th e Chandra error circle are those detected in 
the band bv lWachter et"ai] fcooa) . 



6 OPTICAL OBSERVATIONS WITH THE GRAN 
TELESCOPIO CANARIAS 

We observed the field of SGR0418-F5729 through the Sloan i'-band 
filter (central wavelength A = 7800 A) on 2009 September 15 (see 
Fig. [TJ with the Optical System f or Imaging and low Resolution 
Integ rated Spectroscopy (OSIRIS: iRasilla. Lopez Lopez & Teiadal 
I2OO8?) camera mounted at the 10.4-m Gran Telescopio Canarias 
(GTC) at the Observatorio del Roque de los Muchachos in the La 
Palma island (Spain). OSIRIS is a two-chip CCD detector covering 
the wavelength range between 3650 A and 10500 A, with a nomi- 
nal fov of 7.8 X 8.5 arcmin^ which is actually decreased to 7 x 7 
arcmin^ due to the vignetting of chip #1. The unbinned pixel size 
of the CCD is 0'.'125. 

We took 25 dithered exposures of ~ 100 s each, for a total inte- 
gration time of ^--^2500 s, under dark time and with average seeing 
conditions of ~l'.'3 (air mass 1.3-1.5). SGR 0418-1-5729 was po- 
sitioned in chip #1, slightly offset from the CCD centre to avoid 
occultations from the vignetted regions during the dithered expo- 
sure sequence. To increase the signal-to-noise ratio, the exposures 
were taken with a 2 x 2 binning. Bias and dark frames were taken 
as daily calibrations. In particular, dark frames were taken to cor- 
rect for the anomalously high dark current affecting the CCD ex- 
posures and caused by technical problems. Due to bad weather 
conditions, no sky flats were taken on 2009 September 15 and 
those taken on 2009 September 5 were used instead. Data reduc- 
tion (bias and dark subtraction, flat-field correction) was performed 
using standard tools in the IRAF packag^ll CCDRED after trimming 
each frame for the vignetted regions. Single dithered exposures 
were stacked and averaged using the task DRIZZLE which also 
performs the cosmic ray filtering. Photometry calibration was per- 
formed using observations of the standard star field PG 1528-1-062 

See |http://iraf .noao.edu/ 1 



( lLandoll|[l983h . The astrometry of the science frames was cali- 
brated against the pos itions and coord inates of Two-Micron All- 
Sky Survey (2MASS: Isknitskie et alJi2006 ) stars, resulting in an 
overall accuracy of ^CK'! in both Right Ascension and Declina- 
tion. 

The computed SGR 04 18-1-5729 position on the OSIRIS im- 
age is shown in Fig. (6) We did not detect any source within 
the 95%-confide nce Chandra erro r circle (radius of 0'.'35) of 
SGR 04 18-1-5729 jWoods et alj2009l) . In particular, we did not find 
evidence for the possible detection of the infrared source (Ks — 
21.6 ± 1.3 mag) tentatively seen by Wachter et al. (2009) south- 
west of the Chandra error circle on 2009 August 2, using the Wide- 
field InfraRed Camera (WIRC) of the Hale 5.1-m telescope at Palo- 
mar Observatory (California, USA). We set a 3(j limiting magni- 
tude ^f_i^_^_^25J_on the optical brightness of SGR 0418-1-5729 
(Mig nanietaLll2009l) . which is the deepest obtained so far in the 
optical band for this source. The previous optical upper limit, ob- 
tained from observations performed with the Auxiliary-port Cam- 
era (ACAM) Imager at the 4.2-m William Herschel Telescope 
(WHT) on La Palma o ne month earlier, on 2009 August 16, was 
in fact r > 24 mag jRatti et ai]|2009l) . After de-reddening our 
measured upper limit using the value of the hydrogen column den- 
sity derived from the X-ray spectral fit|f| (Table [2j, and using the 
unabsorbed 1-10 keV X-ray flux measured by Swift/XRT (obser- 
vations 008-009, which were carried out at an epoch close to 
the GTC observations), we infer an X-ray-to-optical flux ratio of 
F-x/Fi, > 10 000 (for a bandwidth A A = = 1500 A), in l ine with 
the ratio observed for other magnetars (e.g. lMignanill2009l) . 

Similarly, we constrained the X-ray-to-near-infrared flux ra- 
tio using as a reference the Hale telescope observations, for which 

^ We u sed the relation between the A^h and the interstellar extinction Ay 
of iPred ehl & Schmitt ( 1995), a s well as the relati ons between the extinc- 
tions at different wavelengths of lFitzpatrickl il999l) . 
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we assumed a non-detection of SGR 0418+5729 down to Ks = 
22.9 mag. Since no Swift/XKT observations were obtained quasi- 
simultaneously to these observations, we assumed the 1-10 keV 
X-ray flux obtained from the monitoring RXTEfPCA observations. 
After de-reddening, we inferred Fx/Fk^ > 50000. The X-ray- 
to-optical flux ratio upper limit derived from WHT observations 
is Fx/Fr > 3000. The uncertainties on these ratios, mainly due 
to the uncertainty in the RXTEfPCA flux determination, are of the 
order of Ri20%. 



7 PRE-OUTBURST X-RAY OBSERVATIONS 

There are two Rontgen Satellite (ROSAT) All-Sky Sur- 
vey jVoges et alj ll99S ) datasets containing the region of 
SGR0418-I-5729 (sequences rs930707n00 and rs930708n00); 
the data were collected in the period between 1990 August 8 and 
1991 February 18. The field of SGR 0418-1-5729 received each time 
roughly 530 s o f exposure with the Posit ion Sensitive Proportional 
Counter (PSPC; |Pfeffermann et al J 1987b . 

SGR 04 18-1-5729 was not detected in these data and the upper 
limit on the count rate, computed as the 3a noise level due to the 
background at the SGR location, is about 0.05 counts s~^ in the 
0.4-2.4 keV energy range. Using the WEBPIMM^ (v3.9j) tool 
and assuming an absorbed blackbody with kT = 0.8 keV and 
A''h ~ 1.1 X 10^^ cm^^ (see TableO, the limit translates into an 
absorbed 1-10 keV (2-10 keV) flux of Ri3 x 10"^^ erg cm~^ s~^ 
(«2 X 10"^^ ergcm-2 s"^). 



8 DISCUSSION AND CONCLUSIONS 

A short period of bursting activity recently led to the discovery of 
the soft gamma-ray repeater SGR 04 18-1-5729 Ivan der Horst et all 
l2009l :l van der Horst e t al. 2010). In the first days after the onset of 
the outburst, the observed flux was of a few 10"^^ erg cm^^ s^^, 
and it decreased by a factor ~10 during the following five 
months. It is unknown whether the flux of SGR 04 18-1-5729 has 
already decayed to (or dropped below) the pre-outburst level, 
since the archival (ROSAT) data of the field provide only an up- 
per limit of ~2 X 10^^^ erg cm^"^ s^^ and magnetars often dis- 
play flux variations over an order o f magnitud e or more (see e .g. 
iGotthelf & Halpern 2007: llsrael eta l. 2007; Es posito et all2008l) . 

The X-ray data presented here allowed us to measure with 
accuracy the spin period of SGR 04 18-1-5729 and to reassess the 
results reported in llsrael et alj ( |2009|) . where a (marginally signif- 
icant) period derivative of P = 3.0(9) x 10"^'^ s s^^ was pro- 
posed. Now, on a longer baseline, we only infer a phase-coherent 
3a limit on the period derivative of |P| < 1.1 x 10~^^ s s~^ (valid 
ove r the range MJD 54993-55144; see iKuiper & Hermsed [2009 
and lWoods. Gogus & Kouveliotoull 20091 for previous upper limits). 
The corresponding upper limits on the surface dipole magnetic field 
stre ngth (inferred within the usual vacuum dipole framework; see 
e.g. iLorimer & Krameij|2004h and spin-down luminosity are B ~ 
3.2xlO^®(PP)^/^ < 3x10^^ G and £; = 47r2/PP-=' < 6x10^" 
ergs~^, respectively, where I ~ 10*'' g cm^ is the moment of 
inertia of a neutron star. The properties inferred from the period 
derivative should be taken with additional caution for magnetars, 
since their spin-down rates can be highly variable. Nevertheless, 

^ See |http://heasarc.gsfc.nasa.govArools/w3pimms.html| 



the present limits on B and E for SGR0418-H572 9 are the low- 
est va lues among magnetars (see e.g. Figure 3 of lEsposito et al.l 
l2009h . IE 2259-1-586 bein g the cl osest, with P ~ 6 x 10^^ G and 
P ~ 5.6 X 10^^ erg s'^ [Gavriil & Kaspi 2002). 

The distance of SGR 0418-1-5729 is not known yet; how- 
ever, considering the structure of the Galaxy proposed by 
IHou. Han &"shil l l2009l) . the Galactic coordinates of the SGR (/ = 
147? 98 and b — 5? 1 2) suggest that it is locate d within a distance 
smaller than ~5 kpc. Ivan der Horst et al J f2O10l) assume a distance 
of ~2 kpc, in the hypothesis that SGR 0418-1-5729 is situated in the 
Perseus arm of the Galaxy. Such a relatively small distance is con- 
sistent with the measure of the photoabsorption derived from the 
X-ray spectral analysis (TVh ~ 1.1 x 10^^ cm^'^, see Table[2}. This 
value is in fact significantly smaller than the total Galactic interstel- 
lar hydrogen alone the line of sight which can be estimated from the 
Leiden/ Argentine/Bonn (LAB) survey of HI ('Kalberla et al .ll2005h 
and the CO survey of Dame, H<u;tmann & Thaddeus ( 200 1 ) for the 
value of molecular H2 as A^h = A^hi + 2iVH2 = (4.30 + 2 x 
0.67) X 10^^ cm"^ ~ 5.6 x 10^^ cm'^ 

Conservatively assuming a distance of 5 kpc, the 1-10 keV 
X-ray luminosity measured with Swift varied during the outburst 
from Lx ~ 4 x 10^* erg s^^ to 6 x 10^^ erg s^^. The minimum 
observed ratio between the luminosity and the spin-down energy 
loss is Lx/E > 10^d|, where dN indicates the distance in units of 
N kpc, a value which, when compared to those of other magnetar 
sources, fully qualifies SGR 0418-1-5729 as a member of the mag- 
netar class. 

The upper limit on the period derivative corresponds, in the 
standard magneto-dipole braking model, to a characteristic age 
Tc = P/ (2P) > 1.1 million years, the highest (so far) among mag- 
netars. The location of SGR 0418-1-5729 on the P-P diagram (see 
e.g. Figure 4 of .McLaughlin et al. 2009 for an updated diagram) 
falls quite close to the region of the X-ray dim isolated neutron 
stars (XDINSs; see lTurollall2009l for a review), often proposed to 
be evolved magnetars on the basis of their comparable spin periods 
and relatively large dipole magnetic fields. Although age estimates 
based on the magneto-dipole model are highly uncertairQ it is in- 
triguing to think SGR 04 18-1-5729 as an 'old magnetar' (or a young 
XDINS), ideally linking the two classes of isolated neutron stars. 

We performed the spectral analysis of the Swift observations 
of SGR0418-F5729 by assuming different phenomenological mod- 
els (combinations of blackbodies and power laws) or more com- 
plex ones, based on resonant cyclotron scattering of seed surface 
photons by magnetospheric currents (ID and NTZ models). In all 
cases, we found that Swift data are compatible with a scenario in 
which the flux decrease is mainly due to the cooling and shrinking 
of a hot emitting region on the neutron-star surface, which may 
have been suddenly heated following the bursting activity. This 
is in agreement with what is observed in the post-outburst evolu- 
tion of other transient magnetars (e.g. Gotthelf & Halpern 2005]). 
SGR 0501-1-4516 is the only other repeater for which the outburst 
decay has been monitored with detailed, quasi-continuous observa- 
tions (XMM-Newton and S wift) over a time span (^^^160 days) com- 
parable to the present one (Rea etalj|200^ . In that case, however, 
the flux decline was exponential with a timescale ^^24 days and 
no evidence for a power-law decay. The first Swift observation was 
performed only 33 days after the onset of the bursting phase and 



An independent age estimate could come from the association with a 
supernova remnant. However, there are no known remnants within 2° from 
SGR 0418-1-5729. 
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by that time one may expect that majo r deviations from a d ipolar 
magnetosphere have died away (see e.g. lEsposito et alJlOOSi) . leav- 
ing only a cooling hot spot at the neutron-star surface as the relic 
of the outburst. Alternatively, the twisted portion of the magneto- 
sphere could have been quite small, as suggested bv lSeloborodovl 
( l2009h in the case of the AXP XTE J1810-197, making it difficult 
for thermal photons to resonantly scatter onto magnetospheric cur- 
rents. 

RXTE pulse profiles (and, to a lesser extent because of the 
lower count statistics, also Swift ones) exhibit a complex pattern, 
with an overall double-peaked shape. If we assume that emis- 
sion mainly comes from the star surface with little contribution 
from the magnetosphere (see the discussion above) and it is es- 
sentially isotropic|j double-peaked light curves arise when the 
angles 5, between the line connecting the two (antipodal) emit- 
ting spots and the rotation axis, and x> between the line-of-sight 
and the rotation axis, are such that £, + x 120° for a ra- 
dius Ans ~ ?>Rs (here Rs is the Schwa rzschild radius; e.g. 
lBeloborodovll2002l J Poutanen & Beloborodovll2006) . This suggests 
that SGR 0418-1-5729 is a nearly orthogonal rotator seen at a 
large inclination angle. Pulse profiles are both energy- and time- 
dependent, with evidence for a hard (>4 keV) component con- 
centrated near the second light maximum (phase ~0.6), which be- 
comes less prominent as the flux declines (see Fig. [3}. This, again, 
points to a picture in which only a limited portion of the star surface 
has been heated as a consequence of the event which produced the 
outburst. The near coincidence in phase of the (second) soft and 
hard peaks, together with the decrease of the latter at later times, 
may indicate that heating involved only one of the two spots, e.g. 
producing a hotter region inside (or close to) one of the warm polar 
caps. 

Phase-resolved spectroscopy of the (combined) Swift WT data 
shows that emission near the phase intervals of maximum soft- 
ness/hardness is described well by the superposition of two black- 
bodies with different radii and temperatures in each interval. The 
radii of the two cooler blackbodies are not much different and, al- 
though somewhat too small (~ 4-6 km), might be associated with 
the neutron-star radius, while their temperatures are quite close 
and consistent with being the same within the errors. A possibil- 
ity, then, is that emission comes from two spots at different tem- 
perature, while the rest of the surface is at a lower temperature. 
Phase-resolved spectral fits suggest two other options: two emit- 
ting regions, each comprising a hotter inner cap and a surrounding 
warmer corona, and two equal spots, one of which emits a black- 
body plus a Gaussian absorption feature. In both these scenarios 
the rest of the surface is likely to be rather colder than the lower 
temperature obtained from the blackbody fit. A preliminary cal- 
culation indicates that all these configurations can produce pulse 
profiles which are in general agreement with the observed ones. 

After the outburst, the flux of SGR 04 18-1-5729 follows a bro- 
ken power-law decay, F oc t", with a break at ~19 days and 
ai ^ —0.3, 02 ~ —1.2 (for the fit to the combined RXTE and 
Swift data; see Section 13. 2t . The post-outburst cooling of mag- 
netars has been modelled in terms of a sudden heat deposition 
in the neutron- star crustal layers, due to the release of energy 
stored in the toroidal component of t he inter nal magnetic field 
jLvubarskv, Eichler & ThomDsonl2003 ; see also lKouveliotou et al.l 



° We warn that this is an oversimplification since thermal emission from a 
magnetised neutron star, either from an atmosphere or a cond ensed surface, 
is not isotropic fe.g. lZavlin et alJl995l ; [Turolla. Zane & Drakei2004l) . 



l2003h . In this picture the flux decays as F oc t~"^'\ with n ~ 2- 
3. While the predicted index might account for the flux evolu- 
tion in SGR 0418-1-5729 after the break, it is inconsistent with 
what is observed before the break, and, even more important, the 
model does not predict any change of slope during the decay. 
On the other hand, it should be stress ed that the calculation of 
iLvubarskv. Eichler & Thompson! ( l2002h was mainly aimed at ex- 
plaining the flux decay in SGR 1900-1-14 after its giant flare and 
consequently assumes a magnetic field B ~ 10^^ G. According 
to the limit on P presented in Section [4] the magnetic field of 
SGR 04 18-1-5729 is likely to be two orders of magnitude lower, so 
the model may not be directly applicable in this case. 



A different possibility is that the surface layers of the star are 
heated by returning curren ts that flow along the closed field lines of 



a twisted magnetosphere i 


Thompson. Lvutikov & Kulkamill2002l; 


iBeloborodov & Thompson 


2007; Beloborodov 20091). The twist is 



likely to be localised into a small bundle of current-carrying field 
lines (the j-bundle) close to the magnetic pole and must decay in 
order to supply its own supporting currents. As the magnetosphere 
untwists, the angular extent of the j-bundle decreases , the luminos- 
ity decays and the area hit by the currents shrinks j Beloborodovl 
2009). The details of the evolution depend on the initial distribu- 
tion of the twist inside the j-bundle and no general law for the 
luminosity decay exists. If the twisted region is small, the emer- 
gent spectrum is dominated by emission from the heated cap and 
it is thermal because photons do not have many chances to scat- 
ter. The flux is then F oc R^T*^ where R and T are the radius 
and temperature of the cap. The decay of the temperature (as de- 
rived from the single blackbody fits to RXTE and Swift data) is 
reasonably well described by a single power- law, T oc t"" ''^ (see 
Fig. [TJ, over the entire time span, without the need to introduce 
a break. If the temperature variation is entirely due to the cool- 
ing of a cap, the emitting area decreases as oc f^"-''** and 
R? oc f""'^^ before and after the break, respectively. If the sin- 
gle blackbody model mostly traces the emission from the hot re- 
gion, as it seems reasonable also from the comparison of the fit pa- 
rameters with those of the hotter component in the two blackbody 
model (see Tableland Section[5j, it follows that, after the break, 
the radius of the cap scales as f^"'^^, close enough to the depen- 
dence derived above (R oc t^'^'^*). A better estimate of the cap 
radius could come from the two blackbody model, but the limited 
statistics allows its application only to the combined Swift obser- 
vations, taken at different epochs after the break. Assuming as a 
rough estimate that the hottest blackbody radius R2 ~ 0.43 km 
can be associated with t ~ 35 days, around which five observa- 
tions cluster, the source luminosity is L ~ 4 x 10'^'*d| ergs~^. 
This value is too high to be explained by the rate of Ohmic dissi- 
pation, L « 10^''(g/10^^ G)(i?Ns/10'^ CYa)iP{V/l(f V) sin* 9, 
erg s "^teeloborodovl'2009'; here i> < 1 rad is the twist angle and 
V is the discharge voltage) if the magnetic field is ~ lO'^^ G and 
the angular size of the twist is sin 6'* ~ R2/RNS ^ 0.03. The 
problem is however eased for a larger cap area, as in the case in 
which the cap contains a small hot region surrounded by a warm 
corona. 

Currently, several AXPs and SGRs have be en firmly identi fied 
also at optical and/or infrared wavelengths (see lMignanill2009l and 
references therein). However, there are no detailed predictions yet 
about the X-ray-to-optical (or infrared) flux ratio that a magnetar is 
expected to have. Our deep GTC i' observation of SGR 0418-1-5729 
points to an X-ray-to-optical flux r atio exceeding 10"*, a value not 
unheard of among magnetars (e.g. iMignanillioogI) . The GTC im- 
age adds a new piece to the still sparse - but growing - database 
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of optical/infrared observations of magnetars exploring a range of 
different states of high energy activity. Enlarging such a database is 
indeed a farsighted effort since it will be the only way to decipher 
the nature of magnetar optical/infrared emission and of its correla- 
tion with the high energy emission. 
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